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DERIVATION OF APPROXIMATE EQUATIONS FOR SOLVING 
THE PLANAR RENDEZVOUS PROBLEM 

By Gene W. Sparrow and Douglas B. Price 
Langley Research Center 

SUMMARY 

The desirability of obtaining simplified and accurate relative equations of motion 
for two-impulse rendezvous application resulted in the examination of the Clohessy and 
Wiltshire type of simplified equation. An analysis of these equations revealed that large 
errors were possible if the radial separation of the two vehicles was large. In an effort 
to develop more accurate equations for rendezvous, the apparent erroneous terms were 
replaced by exact terms. The resulting modified equations contained terms not present 
in the Clohessy and Wiltshire equations. When the modified equations were applied to 
the problem of rendezvous, they were found to be generally superior to the Clohessy and 
Wiltshire equations, particularly for long transfer angles. The scope of the investigation 
is limited to the case in which the target vehicle is in a circular orbit coplanar with the 
ferry vehicle's orbit. 


INTRODUCTION 

Of particular interest in the space program are the problems of rendezvous 
between two vehicles in orbit. Since the exact equations describing the relative motion 
of the vehicles are nonlinear, the usual approach in finding a solution is to approximate 
the exact equations. Clohessy and Wiltshire first obtained a linear approximation of 
these equations for the case in which the target vehicle is in a circular orbit (refs. 1, 2, 
and 3). Some forms of the linear equations that are even simpler, though less accurate, 
were obtained by modifying the gravity term (ref. 3). London (ref. 4) improved the accu- 
racy of the linear form by including second-order terms in his equations. Anthony and 
Sasaki (ref. 5) and others generalized further by considering a slightly elliptical target 
orbit. 

To obtain the impulses required for a two-impulse rendezvous, it is necessary to 
know the initial required velocity and the velocity at rendezvous. The linear equations of 
Clohessy and Wiltshire are sufficiently simple so that the two impulses may be defined 
explicitly in closed form. Although it is possible to define explicitly the two impulses 
for rendezvous in a second-order system (see ref. 5), the expressions are usually not 



simple. Since the linear form is desirable for simplicity but is not accurate, an effort 
was made to find more accurate equations of this form for two-impulse rendezvous 
application. This paper presents the results of the investigation. 

The scope of the investigation is limited to the case in which the target vehicle is 
in a circular orbit coplanar with the ferry vehicle’s orbit. The first impulse is assumed 
to occur at the initiation of the problem, and the target vehicle's angular travel during 
the coast following the first impulse is arbitrarily limited to 360°. Examples are pre- 
sented which illustrate the application of the equations derived herein to the problem of 
rendezvous. 


SYMBOLS 


d 

d R 

k 

Q(y 0 ) 

r F 

r T 

t 


slant-range separation of ferry and target vehicles, kilometers 
slant-range separation of ferry and target vehicles at intercept, kilometers 
an arbitrary integer 
function defined by equation (13) 

ferry vehicle's distance from mass center of central body, kilometers 

target vehicle's distance from mass center of central body, kilometers 

elapsed time from problem initiation, seconds 

shell-coordinate reference-axis system (see fig. 1) 

coordinates in the X,Y axis system, kilometers 

inertial reference -axis system fixed in central body (see fig. 1) 

error expressions defined by equations (12) 

magnitude of first velocity impulse, kilometers/second 
magnitude of second velocity impulse, kilometers/second 
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6 symbol used to denote a trigonometric expression defined by equation (17) 

l± product of universal gravitational constant and mass of central body 

Up angular position of ferry vehicle as measured in Xj,Yj axis system, 

radians 

r elapsed time from first impulse to intercept, seconds 

4> angle defining the direction of the ferry vehicle with respect to the target 

vehicle, degrees (see fig. 1) 

cu angular rate of target vehicle, radians/second 

cut angular position of target vehicle as measured in the Xi,Yi axis system, 

radians or degrees (see fig. 1) 

cur angular travel of target vehicle about central body during the interval 

extending from first impulse to intercept, radians or degrees 


Subscripts: 

c circular component 

o initial conditions 

An asterisk denotes a required initial condition. 

Dots above symbols denote differentiation with respect to time or with respect to 
cut when nondimensional form is used. 

A bar above a symbol denotes the nondimensional form. 

ANALYSIS 

Coordinate System and Equations of Motion 

The equations of relative motion between the ferry and target vehicles are to be 
expressed in so-called shell coordinates. (See fig. 1.) The shell coordinate system was 
first introduced in reference 6. In this system the X and Y axes are centered in the 
target vehicle, which is moving in a circular orbit at a radius rrp with angular rate cu. 
The ferry vehicle's position relative to the target vehicle is defined by the x-coordinate, 
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which is measured along a curved line of constant radius rp, and by the y-coordinate, 
which is directed along r F and measured from the target vehicle's radius r T . The 
rendezvous maneuver is initiated at first impulse (cot = 0) with the target vehicle located 
on the Xi axis and is completed when x = y = 0 for some arbitrary value of cot. The 
symbol d denotes the slant-range separation of the ferry and target vehicles at any 
time. 

The inertial positions of the ferry and target vehicles are defined by the polar 
coordinates r F ,Op and r-p, cot, respectively. The relationships between these polar 
coordinates and the shell coordinates are 

x = r-p^cot - a F ' 
x = r T (co - a F ) 
x = -rpdp 



y = r 


( 1 ) 


The differential equations of motion of the ferry vehicle about a central body in polar 
coordinates are expressed by the familiar form 




r _ r a 2 _ _ M 

r F r F a F “ 2 


r F^F + ^F°F ” 0 


( 2 ) 


The differential equations of motion in shell coordinates may be determined by substi- 
tuting into equations (2) the appropriate terms found from equations (1). The exact dif- 
ferential equations of motion in shell coordinates thus obtained are 



( 3 ) 


The gravitational constant | u has been eliminated from equations (3) by the substitution 
ju = rrp3a>2 ) s i nC e rpco = r-p for a circular target orbit at rp. 

It is convenient to use equations (3) in a nondimensional form which is obtained by 
nondimensionalizing displacements with respect to r-p, velocities with respect to the 
target velocity corp, and accelerations with respect to co^rp. The nondimensional 
form of equations (3) is 
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x = 2(1 + y) -1 (l - x)^ 
y = (1 + y)(l - x) 2 - (1 + y) 


( 4 ) 


Approximate Differential Equations of Motion and Their Solution 
1 2 / • \2 

When the terms (1 + y)~ , (1 + y)“ , and (1 - x) in equations (4) are expanded 
in a binomial series and grouped according to their order, the following equations are 
obtained: 


x = 2yjj. - (x + y) + (yx + y 2 ) - (y 2 x + y 3 ) 
y = (-2x + 3y) + (x 2 - 2yx - 3y 2 ) + (yx 2 + 



(5) 


If it is assumed that the ferry vehicle remains in the vicinity of the target vehicle during 
the rendezvous maneuver, the higher order terms of equations (5) will be small in com- 
parison with the lower order terms and may be neglected. A first-order approximation 
of equations (5) is 



( 6 ) 


Equations (6) are the so-called nondimensionalized Clohessy and Wiltshire equations in 
shell coordinates. Equations (6) may be easily integrated with respect to time to obtain 


x = (4 sin cot - 3cot)x 0 + 2(1 - cos cot)y Q + 6(cot - sin u>t)y 0 + x 0 ( 
y = -2(1 - cos cot)x 0 + (sin cot)y 0 + (4 - 3 cos cot)y 0 


(7) 


where the initial conditions x 0 , y Q , x 0 , and y 0 are evaluated at cot = 0. These 
equations are a linear first-order solution to the exact equations of motion. Equa- 
tions (7) are readily applied to a two-impulse rendezvous problem since it is possible to 
define explicitly the required initial velocity for intercept and the resulting velocity at 
intercept. However, equations (7) are valid only for small values of |y| and small 
transfer angles cot. 

Since equations (7) permit closed-form solutions, it was felt that an effort must be 
made to improve this type of equation for rendezvous application. To this end, equa- 
tions (7) were examined in a special way in an effort to define the cause of their inaccu- 
racy. A possible way to improve the equations became evident after examination. 
Although the improvement lacked mathematical rigor, the approximate equations 
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(obtained from eqs. (7)) were found to be generally more accurate than equations (7) for 
rendezvous application. 

Examination of the Approximate -Solution Equations 

The rendezvous trajectory may be considered as a trajectory which has been per- 
turbed about a circular orbit at the initial relative altitude y = y Q . If the eccentricity 
of the trajectory is small, the circular part is of primary significance, and it is impor- 
tant to know how this part is mathematically expressed by equations (7). This may be 
determined by first considering x 0 to be composed of two parts, with one part x c 
being the value of x 0 for a circular orbit at y 0 . Thus 

x 0 = x c + Ax 0 (8) 

Substituting equation (8) 
lowing equations: 

x = 


y 


The terms are arranged so that if the initial conditions correspond to a circular orbit at 

R ♦ 

the initial altitude, which requires that Ax 0 = y Q = 0, the second bracketed terms of the 
x and y equations vanish. The first bracketed terms then determine in a linear sense 
x and y as a function of cot for the ferry being initially in a circular orbit. The 
terms are functions of x Q and y Q (both of which may presumably be measured), x c 
(which is computed), and cot. 

Exact expressions for x c , x c , and y c under the circular orbit assumptions are 

x r = x r cot + x n l 

> <10a) 

y c = y 0 J 

• / . v-3/2 

x c = 1 - [1 + y Q J (10b) 


into equations (7) and rearranging the terms leads to the fol- 


|(4 sin cot - 3cot)x c + 6(cot - sin cot)y Q + xjj 
+ J^4 sin cot - 3cot)Ax 0 + 2(1 - cos cot)y^j 


r 




j-2(l - cos cot)x c + (4 - 3 cos cot)y ( 


+ -2(1 - cos cot)Ax Q + (sin cot)y 


?o 


(9) 
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Equation (10b) is determined from equations (4) by setting y equal to 0. When com- 
paring equations (10a) with the first bracketed terms of equations (9) it is difficult to see 
the error source. A more meaningful comparison between these equations is possible 
if x c is expressed as a function of y Q . By expanding (l + in a binomial 

series, equation (10b) may be written as 


x c =fy 0 - Q(y 0 ) (u) 

where Q^y 0 ) * s used to represent the second-order and higher order terms of the 
expansion. 

When equation (11) is substituted into equations (10a) and into the first bracketed 
terms of equations (9), and the latter equations are subtracted from the former, expres- 
sions for the error are obtained: 


(*c) e rror = 4Q W (sin Wt ' Wt) , 

(terror = ’ 2Q (^) H ’ cos "‘i 


( 12 ) 


The error thus depends upon the higher order terms of y Q and of a>t. No error is 
present if the ferry vehicle is initially at the same altitude as the target vehicle. The 

(x c ) term is noncyclic and the error increases with wt. The (y„) term is 

cyclic and the error is zero for transfer angles that are multiples of 2 t \ radians. The 
maximum value of (yc) error occurs at (2k - 1) 7 r radians (k is an integer). Although 

it is possible to approximate Q^o) by using only its lowest order term ^ y Q ^, the 

results presented in this paper are based upon the exact value of Q(y 0 )- The exact 
expression is 


Q(y 0 ) = f y 0 - 1 + ( x + yo ) 3/2 ( 13 ) 

which is determined by equating equation (10b) to equation (11). 

Some example calculations will better illustrate equations (12). A lunar orbit is 
assumed, with the reference or target orbit located at an altitude of 148.16 kilometers 
(80 nautical miles) above the lunar surface. The ferry vehicle is initially located at dif- 
ferent altitudes from the target vehicle and its initial velocity is that required for a cir- 
cular orbit. The errors computed from equations (12) for these conditions are illus- 
trated in figure 2, for six initial altitudes of the ferry vehicle. Although the case 
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y Q = -250 kilometers is physically impossible, it is included for completeness. Fig- 
ure 2(a) shows that (x c ) err0 r increases sharply with both y Q and cot, and that the 

errors are relatively large for the longer transfer angles. Although the (Yc) error 
term, figure 2(b), is smaller than the corresponding (xc) error term, it is still signifi- 
cant. On the basis of equations (12) and figure 2, the predicted values of x and y for 
any value of cot tend to be larger than they should be. For example, the ferry vehicle 
initially located in a circular orbit at y 0 = -100 kilometers will rise to an altitude of 
-60 kilometers at cot = 180°. The value of x is about 135 kilometers greater than it 
should be for this case. 

Figure 2 has shown that large errors occur in the first bracketed terms of equa- 
tions (9) for large transfer angles when a circular orbit is assumed. Equations (9) when 
applied to rendezvous are also inaccurate for large transfer angles. These two facts 
suggest that the inaccuracy of equations (9) may be due to the inaccuracy of their first 
bracketed terms, particularly when y Q is large. This line of reasoning also suggests 
that better results may be obtained for large transfer angles if the inaccurate first 
bracketed terms are replaced by the exact terms of equations (10a). This approach will 
be used in deriving a new set of equations of relative motion. 

Derivation of Modified Equations of Motion 

The first bracketed terms of equations (9) are replaced by their exact expressions 
from equations (10a), and Ax 0 is eliminated from the second bracketed terms of equa- 
tions (9) through the use of equation (8). The resulting equations depend upon the initial 
conditions x 0 , y Q , x 0 , and y 0 , and upon x c . The x c terms may now be replaced 
by equation (11) so that the resulting equations will depend upon only the initial-condition 
terms and cot. When these operations are performed and the equations are simplified, 
the new modified equations of motion are obtained: 


x = (4 sin cot - 3cot)x 0 + 2(1 - cos cot)y 0 + 6(cot - sin cot)y 0 + x 0 + 4Q^y 0 ^(sin cot - cot)| 
y = -2(1 - cos cot)x 0 + (sin cot)y Q + (4 - 3 cos cot)y 0 - 2Q^y 0 j(l - cos cot) 


(14) 


where Q(y 0 ) is defined by equation (13). Equations (14) are the solution to the differ- 
ential equations 


x = 2y 
y = -2x + 



( 15 ) 


The difference between equations (15) and equations (6) may be noted. 
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Equations (14) differ from the original equations (7) by including the terms 
involving Q^y 0 ), which depend only on y Q and cot. 


Derivation of Two-Impulse Equations for Rendezvous 


In order to determine the first and second impulses, it is necessary to know the 
initial required velocity components for intercept and the final velocity components at 
intercept. The initial required velocity components x Q and y Q are determined by 
setting x = y = 0 in equations (14) (intercept condition) and solving for the initial con- 
ditions x Q and y Q . The following equations are obtained: 




sin cor + 6) 


x Q = i^o sin WT + y Q (6coT sin cor -14 + 14 cos coT)j - ^l(wr 
y°* =|[2 (1 - cos cot)x q + y 0 (4 sin cot - 3cot cos cot)J + 0(yo) (6 

- 3cot sin cor) 


> (16) 


where 


6 = 3cot sin cot - 8(1 - cos cot) 


(17) 


Equations (16) determine the initial required velocity components for intercept as a 
function of the initial displacement x 0 ,y 0 anc * the desired transfer angle to intercept 
cot = cot. The velocity components at intercept are determined by using equations (16) 
as the initial velocity components of the trajectory. When equations (16) are substituted 
into the velocity equations defined by differentiating equations (14) with respect to cot, 
the velocity components at intercept may be obtained as 


-i^ si 


x(cot) = - x Q sin cot + 2y 0 (l - cos cot) ^i~(coT sin C0T + 5) 


3 - 


q ( 5 °), 


y (wt) 


* f [- 2S ° 


(1 - cos cot) + y 0 (4 sin cor - 3 cot) j-^-(6 - 3 cot sin cot)cot 


3 - 


46 




(18) 


The first impulse is the vector difference between the initial required velocity for 
intercept and the actual initial velocity. The magnitude of the first impulse is defined 
by the dimensional equation 


AV 1 = 


|[ i ( r F,o i F,o3 2 + [ 4 ( i 'F,o)]‘ 


(19) 
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where A^r^, is the change in the circumferential velocity of the ferry vehicle 

and A/rp A is the change in its radial velocity, or 


A ( r F,o CT F,o) ~ r F,o( ff F,o " CT F,o 

A ( r F,o) = r F,o " r F,o 
Substitutions from equations (1) lead to 




( 20 ) 


A ( r F,o d F,o) " r F,o 

a (^f,o) = V - y Q 


CO - 


CO 


= -( 1+ 4r°*'*o) 


> ( 21 ) 


The first-impulse equation becomes, in nondimensional form, 


AV X = 


(! + y 0 )( 5 o* - xq) + y D * - yo 


2 n * • 2 


( 22 ) 


Equation (22) represents the actual magnitude of the first impulse applied, expressed in 
shell coordinates. Since the relative velocity between the ferry and target vehicles 
must be brought to zero for rendezvous, the second impulse is simply the vector differ- 
ence between the zero velocity desired and the intercept velocity. The magnitude of the 
second impulse in nondimensional form is 


AV 2 


1 p* "1 

2 

r* n 

X(WT) 

+ 

y(wr) 


(23) 


Since the second impulse is executed at y = 0, a term involving y Q does not appear in 
the second-impulse equation as it did in the first. This -statement may be verified by 
the procedure used in determining equation (22). 


RESULTS AND DISCUSSION 


In order to demonstrate the effectiveness of using the new approximate equations, 
numerous intercept trajectories were computed. The ferry vehicle was assumed to be 
initially located at some slant-range distance d Q from the target vehicle. For each 
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given d D value, the ferry vehicle's initial direction with respect to the target vehicle 
(defined by angle <p Q , fig. 1) was varied from 0° to 360° in increments of 30°. From 
each initial relative polar coordinate position (d O ,0 O ) °f the ferry vehicle the initial 
required velocity for intercept was computed by two methods for transfer angles from 
30° to 330°. These computed values of required velocity were then used as the input to 
an exact-equation computational program which determined d R’ the slant-range sep- 
aration of the two vehicles at the intended intercept. Method B uses equations (16). 
Method A uses the Clohessy and Wiltshire equations, equations (16) with Q(y 0 J = 0. The 
resulting dR values for both methods were compared. Rendezvous maneuvers were 
assumed to occur about only one central body, the moon. The target vehicle's circular 
orbit was located at an altitude of 148.16 kilometers (80 nautical miles). 

The results of the investigation (for the two methods) are presented in figures 3 
and 4 as semilog plots. The separation at intercept djj is plotted against the transfer 
angle cot for a range of initial polar angles <p Q . For figure 3, a d Q value of 
100 kilometers was used for cp Q values of 0° to 330°. Figure 4 shows the result for 
d 0 values of 50 kilometers and 250 kilometers for a single <p Q of 90°. In figures 3 
and 4 faired curves were used to connect the data points. 

For y D ~ 0 (figs. 3(a) and 3(g), where cf> 0 = 0° and 180°), both methods result in 
about the same values of dR for any given cot. The magnitude of dR is small com- 
pared with the initial d 0 value. As |y Q | increases toward its maximum (<p 0 = 90° 
and 270°), the dR obtained by the use of method A increases sharply with increasing 
cor. Method B, however, does not exhibit this large increase in dR. For all of fig- 
ure 3, with the exception of <p Q = 0° and 180°, method B results in smaller dR values 
than method A over the cot range of 120° to 300°. In most cases shown, this is true 
for even smaller cot values than 120°. In the region where method A results in 
smaller dR values than method B (small values of cot) the differences in dR are 
small. However, in the region where method B is superior to method A (cot not small) 
the differences in dR values are usually large. The values in the <p Q interval of 0° 
to 180° are approximately the same as the values in the interval from 180° to 360°. 

That is, the curves of the cp Q = 30° case are similar to those of the <p Q = 210° case. 

The two methods were also compared for d Q values from 25 to 750 kilometers, 
but the results are not presented here. The comparison between methods A and B 
remained essentially the same as that presented in figure 3. As might be expected, the 
larger the d Q value used, the larger the resulting dR values. This is illustrated in 
figure 4 for d D values of 50 and 250 kilometers with <p 0 = 90°. The curves of fig- 
ure 4(a) are similar to those of figure 4(b), and also similar to those of figure 3(d), 
where d Q = 100 kilometers and cp 0 = 90°. Both methods deteriorate with increasing 
d Q . When d D is increased by a factor of 5, from 50 to 250 kilometers, the maximum 
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dR of method A increases from 60 to 1200 kilometers, a factor of 20. The near steady- 
state interval of method B increases also by a factor of about 20 from about 3 kilometers 
to about 70 kilometers. 


CONCLUDING REMARKS 

The desirability of obtaining simplified and accurate relative equations of motion 
for two-impulse rendezvous application resulted in the examination of the Clohessy and 
Wiltshire type of simplified equation. An analysis of these equations revealed that large 
errors were possible if the radial separation of the two vehicles was large. In an effort 
to develop more accurate equations for rendezvous, the apparent erroneous terms were 
replaced by exact terms. The resulting modified equations contained terms not present 
in the Clohessy and Wiltshire equations. When the modified equations were applied to 
the problem of rendezvous, they were found to be generally superior to the Clohessy and 
Wiltshire equations, particularly for long transfer angles. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 11, 1968, 

127-51-01-02-23. 
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Figure 3.- Separation of target and ferry vehicles at intercept resulting from the use of Clohessy-Wi Itsh i re equations (method A) and 
modified equations (method B), for transfer angles from 30° to 330° and initial orientation angles % from 0° to 330°. d 0 = 100 k 
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Figure 4.- Separation of target and ferry vehicles at intercept resulting from the use of Clohessy- Wi Itshi re equations (method A) and modified 
equations (method B), for transfer angles from 30° to 330° and initial slant-range separation d 0 of 50 and 250 km. <t> 0 = 90°. 
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